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We present a methodology for the investigation of dislocation energetics in segregated alloys based 
on Monte Carlo simulations which equilibrate the topology and composition of the dislocation core 
and its surroundings. An environment-dependent partitioning of the system total energy into atomic 
contributions allows us to link the atomistic picture to continuum elasticity theory. The method is 
applied to extract core energies and radii of 60° glide dislocations in segregated SiGe alloys which 
are inaccessible by other methods. 
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Dislocations are the most common extended defects in 
crystalline solids. A prototypical example in semiconduc- 
tor physics is offered by the generation of dislocations in 
SiGe alloys grown on Si which relieve the stress induced 
by the 4% lattice mismatch between Si and Ge The 
presence of dislocations in SiGe alloys has a strong im- 
pact on their mechanical and electronic properties. Far 
from the dislocation, the induced strain field depends 
only on macroscopic properties of the material, such as 
its shear modulus. However, near the dislocation line, 
the local properties of the material deviate significantly 
both from those of the ideal solid and the predictions of 
continuum elasticity theory. This "core region" of the 
dislocation, accessible only to atomistic calculations, can 
have different structures depending on the material, the 
applied external strain, and other parameters 

Although the core structures and energetics of disloca- 
tions in elemental crystals have been well studied 0] , 
relatively few investigations have linked the atomistic pic- 
ture of dislocated alloy systems to traditional continuum 
descriptions. Blumcnau et al. used an atomistic descrip- 
tion to calculate core structures and energies in perfectly 
ordered alloys, such as SiC Martinelli et al. per- 
formed simulations on SiGe alloys and found that Si-rich 
and Ge-rich nanowires are formed near the dislocation 
core 6J. However, despite these important contributions, 
a well-defined method for extracting core energies and 
radii from dislocations in segregated alloys has not yet 
been proposed. In addition, the importance of alloy seg- 
regation and its relationship to traditional continuum de- 
scriptions including the concept of core parameters is not 
fully understood. 

In view of the crucial role played by dislocations in 
many important applications for SiGe materials, such as 
bulk alloy layers for quantum wells 0] or nanoislands 
for quantum dots d, 01 j it is highly desirable to under- 
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stand the core energetics and related properties. In the 
present work, we develop a framework to calculate the 
core energies and core radii of segregated alloys using 
atomistic calculations, taking a 60° glide dislocation in 
a Sio.sGco.s alloy as a prototype. This is made possible 
by an environment-dependent partitioning of the total 
energy into atomic contributions, which allows us to lo- 
cally probe the energy at the dislocation core, and link 
the atomistic picture to continuum theories. 

Our approach is based on state-of-the-art continuous- 
space Monte Carlo (MC) simulations, employing the 
multi-component empirical potentials of Tersoff [lo| . 
This method, capable of calculating the free energy of 
the system with great statistical precision, has been ap- 
plied with success in similar contexts [ll|, Three 
types of random moves are involved in the MC algorithm: 
atomic displacements, volume changes, and mutual iden- 
tity exchanges between atoms of different kinds, leading 
to topological and compositional equilibration of the sys- 
tem, as described in Ref. [l^ . The simulation supercell 
consists of 40x10x8 6- atom unit cells of the diamond 
lattice, with lattice vectors having [Oil], [110] and [111] 
directions, including a total of 19200 atoms. We checked 
the convergence of our results with respect to the size of 
the simulation cell by running test simulations with dif- 
ferent cell sizes. Epitaxial SiGe alloys grown on Si(lOO) 
are simulated by constraining them to have the lattice 
constant of Si in the supercell, and then allowing relax- 
ation of the lattice vectors only in the [100] direction. 

The most commonly observed dislocations in the dia- 
mond lattice are those corresponding to the {111} slip 
system; the lowest energy perfect dislocation in tetra- 
hedral semiconductors is the so-called 60° glide disloca- 
tion [s^]. Although such dislocations usually dissociate 
to Shockley partials, especially in diamond, epitaxially 
grown alloys are very likely to hold undissociated 60° 
glide dislocations, due to the strain introduced by the the 
lattice mismatch [l^]- We therefore choose to study 60° 
glide dislocations as characteristic prototypes of disloca- 
tions in SiGe alloys. Two such dislocations with opposite 
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FIG. 2; (color online) Contour plots of the average compo- 
sition (left) and average atomic stresses (right) around a 60° 
glide dislocation in a Sio.sGeo.s alloy at 300 K. The disloca- 
tion line is normal to the plane of the plot. Distances shown 
on the axes of the plots are given in nm. 



FIG. 1: (color online) Equilibrated structure, 2 nm around 
the core of a 60° glide dislocation in a Sio.sGeo.s alloy at 300 
K. The dislocation line is normal to the plane of the plot. 
Sites with c < 0.3, having strong preference for Si, are shown 
in darkest color; sites with c > 0.7, having strong preference 
for Ge, are shown in lightest color. 



Burgers vectors are introduced in the cell 14], so that the 
total displacement field away from the dislocation cores 
is zero. The distance between the dislocation cores is 
about 80 A, ensuring that dislocation motion due to mu- 
tual interaction can be neglected for the temperatures of 
our simulations. The Burgers vectors and the lines of 
the two dislocations lie in the [Oil] and [110] direction, 
respectively. 

We begin by equilibrating the unit cell for dislocated 
bulk SiGe at 300 K. Simulations at higher temperatures 
(900 K) show similar composition profiles. After several 
equilibration steps, including identity exchanges, a long 
run 10^° MC steps) follows, over which the average 
occupancy of each site, c, is calculated. This is defined 
as the fraction of the simulation time (MC steps) when a 
site is occupied by a Ge atom over the whole run. Its val- 
ues range between (site always occupied by Si atoms) 
and 1 (site always occupied by Ge atoms). For the bulk 
Sio.sGeo.s alloy at 300 K, the calculated values of c are 
0.5 ±0.02. However, for the cell containing a dislocation, 
there are sites near the dislocation core with values of c 
very close to 0.0 and 1.0, showing strong preference for 
occupancy by Si or Ge atoms, respectively. This is illus- 
trated in Fig. [1] which shows the equilibrated atomistic 
structure of the dislocation core. It can be observed that 
the alloy segregates into Si-rich and Ge-rich cylindrical 
regions on opposite sides of the dislocation line, which 
can be viewed as self- assembled nanowires in the alloy 



To explain this effect we note the similarity to stress 
induced segregation at the (100) reconstructed surface 
ll| . The significant near-surface stress field, resulting 
from surface dimerization, drives smaller atomic volume 
Si atoms to occupy sites under compression, while the Ge 
atoms prefer sites under tension. In the present case, this 
correlation between local stress and composition can be 
probed via the atomic stress, defined as cr^ = —dEi/dVi, 
where Ei is the e nerg y of atom i and Vi is the volume 
available to atom i llj . In order to link the atomistic pic- 
ture with macroscopic theories of elasticity, we define a 
continuous stress field (T(r) = (7^. The sum runs over 
all atoms i that have Cartesian positions such that 
jr. -r| < 3.0 A. In the same manner, we define a contin- 
uous composition field. Fig. [2] shows contour plots of the 
composition and stress fields on a plane perpendicular to 
the dislocation line. There is a clear correlation between 
the two, as Ge atoms prefer to occupy sites under ten- 
sion while Si atoms prefer sites under compression. The 
observed composition contours are very similar to those 
found by lattice Monte Carlo simulations of dislocations 
in alloys a{r) has the same qualitative features as 
the hydrostatic component of the stress tensor around a 
dislocation, a{x,y) = {a^x + (jyy + o-^z)/3 ~ - jifp-, as 
calculated in linear elasticity theory [3]. Similar segre- 
gation results are found for both relaxed and epitaxial 
(lOO)-constrained alloys. 

The significant segregation evident in Fig. [T]will clearly 
infiuence the dislocation core energy and elastic proper- 
ties. It is therefore important to develop a framework for 
defining dislocation core parameters in the presence of 
segregation. Our approach is based on a partitioning of 
the dislocation formation energy into 'elastic' and 'segre- 
gation' components. Consider first the radial dislocation 
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FIG. 3: (color online) Probability density of atomic energies 
for Si (squares) and Ge (circles) in a Sio.sGeo.s alloy at 300 
K. 

formation energy, 

i 

where the summation is taken over all atoms i contained 
in a cyUnder of radius R and length L around the dis- 
location line. Ehuik is the average atomic energy in 
the random alloy and Ei is the energy of atom i. We 
now decompose the formation energy into two terms, 
Ef{R) — Eei{R) + Eseg{R), where the elastic and the 
segregation components are defined as 

i 

S,es(i?) = -^^(S(C„C',) -Swfc). (3) 

i 

Here, we define E{ci,c'j^) as the energy of a site with av- 
erage occupancy Ci while the average occupancy of the 
neighboring sites is c ■ . Eei (R) represents the energy cost 
of including a dislocation in a segregated alloy (i.e. an al- 
loy with a compositional fingerprint which is identical to 
that produced by a dislocation) . Egeg (i?) is the energy of 
segregation per unit length of dislocation (i.e. the energy 
cost of creating the compositional fingerprint from a ran- 
dom alloy). For elemental solids or ordered heteropolar 
alloys, this is either zero or assumes one of two distinct 
values and can thus be taken into account implicitly [5|. 
For a segregated alloy, however, this decomposition is a 
necessary step in linking atomistic and continuum con- 
cepts. 

To calculate the elastic and segregation energies de- 
fined by Eqs. ^ and ([3]), we first evaluate the func- 
tion E{c, c'). This is achieved by equilibrating a random 



Sio.sGeo.s alloy at 300 K using a 4096-atom cubic super- 
cell. All atoms are fourfold coordinated, and the average 
cohesive energy was found to be —4.19 eV/atom. In Fig. 
m we plot the probability distributions P{E) of atom en- 
ergies. For each type of atom (Si or Ge), we observe five 
distinct and narrow peaks. By inspecting the environ- 
ment of atoms whose atomic energies are close to these 
peak values, we conclude that each peak corresponds to 
a different chemical environment. For Si, the first peak 
is at —4.59 eV and corresponds to Si atoms neighboring 
with four Si atoms; this number coincides with the cal- 
culated cohesive energy of Si at 300 K. The second peak, 
at —4.44 eV, corresponds to Si atoms neighboring with 
three Si atoms and one Ge atom. The third peak, at 
—4.29 eV, corresponds to Si atoms neighboring with two 
Si atoms and two Ge atom. The last two peaks corre- 
spond to Si atoms neighboring three Ge atoms (energy 
—4.14 eV) and Si atoms neighboring with four Ge atoms 
(energy —4.00 eV). The five peaks in the distribution 
of atom energies for Ge atoms can be attributed to the 
same bonding environments. The peak energies and the 
corresponding composition of the neighbors are: —4.39 
eV for four Si neighbors, —4.25 eV for three Si and one 
Ge neighbor, —4.11 eV for two Si and two Ge neighbors, 
—3.95 eV for three Si and one Ge neighbor and -3.81 eV 
for four Ge neighbors. The last number coincides with 
the calculated cohesive energy of Ge at 300 K. This anal- 
ysis yields values of the function E{c, c') for c =0 or 1 
and c' = 0.00, 0.25, 0.50, 0.75 and 1.00; we use linear 
interpolation to calculate values of E{c, c') for all other 
values of c and c'. 

We can now evaluate Eei (R) and Egeg (R) as a function 
of R. The strain field contributing to the elastic energy 
term, E^i (R) , will involve the superposition of the strain 
fields from the dislocation and the compositional segre- 
gation evident in Fig. [2l This compositional segregation, 
consisting of spatially separated regions of compression 
and tension, is clearly associated with a strain field of 
a dipole-like character. It will consequently fall off rela- 
tively quickly compared with the long-range dislocation 
strain field. With this in mind, it is appropriate to de- 
fine several spatial regimes associated with the segregated 
dislocation: The first regime is associated with the dislo- 
cation core radius, i?c, where continuum elasticity breaks 
down. For segregated systems an extended core region, 
also exists, where significant segregation is observed. A 
third intermediate regime corresponds to a superposition 
of the compositional dipole and dislocation strains. Fi- 
nally, the far-field region is dominated by the strain field 
for the dislocation alone. With the intention of bringing 
together continuum and atomistic descriptions we there- 
fore propose a parameterization of the form 

EeiiR) ^b^k\n{R/R,)+E,, (4) 

where Rc and Ec are the respective core radius and en- 
ergy associated with a dislocation of Burgers vector b. 



and k should depend on the region of study. In the 
case of an ideal continuum model of an elemental crys- 
tal, compound or ordered alloy, the factor k is formally 
related to elastic constants of the crystal and geometry 
of the dislocation [2?|. In the case of segregated systems, 
k is expected to assume the random alloy value in the 
far-field region. In the intermediate and extended core 
regimes, the dipole segregated region around the disloca- 
tion will modify the energetics and k must be considered 
as an effective parameter which must be determined from 
atomistic calculations. Since the size of our simulation 
cell seems to be limited to the intermediate regime, the 
values of k we calculate should be smaller than the value 
of k for the non-segregated alloy. 

In Fig. m we plot Eei{R), as calculated using Eq. ([2|) 
as a function of Ini? in the intermediate regime. We find 
that the thus defined Eei{R) is indeed a smooth, lin- 
ear function of Ini?, consistent with the parameterization 
scheme of Eq. (|4]) which effectively takes into account the 
influence of the extended core. We can then determine 
the core radius as the distance at which the linearity of 
Eei{R) breaks down. In Fig. 21 for example, this oc- 
curs at ln(i?/A)=1.5. The core energy is then defined 
as Ec — Eei{Rc), while k is determined by the slope of 
the linear region, in accordance with Eq. The error 
in the somewhat subjective determination of ln(i?c/A) is 
about 0.05, giving an error bar of 0.05 eV/A for E^. By 
taking the average over many runs, we determine dislo- 
cation core parameters for the alloy to be fc = 4.5 GPa, 
Rc= 4.5 A and Ec = £;e/(i?c)=0.59 eV/A. As the core 
parameters describe local deformations around the dislo- 
cation line, they are not very sensitive to a macroscopic 
strain field: For the (lOO)-constrained alloy, we find a 
marginaly lower core energy (0.57 eV/A) and core radius 
(4.2 A). Temperature also seems to play a minor role in 
the core energies: at 900 K, the core energy is only 0.03 
eV/A higher than that at 300 K. 

We have evaluated the radial dependence of Eseg{R) 
directly from Eq. ^ and find a relatively weak depen- 
dence on R which can be understood from Figs. [5] and [31 
Terms in Eq. ([3]) will be of opposite sign depending on 
whether E{ci, c'-) is associated with Ge-rich or Si-rich re- 
gions. Consequently, they will tend to cancel in the sum- 
mation over a cylindrical region about the dislocation 
line, reflecting the symmetry of the segregated regions 
(Fig. [2]). Therefore, we find that Eseg{R) can be approx- 
imated by a constant value of the order 0.2 eV/A. Adding 
this value to Eq. ^ shows that the dislocation forma- 
tion energy in segregated alloys can be approximated by 
the continuum energy expression for a dislocation in an 
elemental crystal with a judicious choice of parameters 
appropriate to specific spatial regimes of the segregated 
dislocation. This surprising result indicates that simple 
continuum concepts, including core energies and radii, 
are still of value in situations where alloy segregation 
is significant and that the effects of the compositional 
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FIG. 4: (color online) Elastic energy, E{R), contained in a 
cylinder of radius R and unit length around the dislocation 
core versus the logarithm of R for a 60° glide dislocation in a 
Sio.sGeo.s alloy at 300 K. 

dipole strain can be included specifically in fc. 

We now apply the above framework for determining 
core parameters in segregated systems to provide addi- 
tional insight into segregation energetics by direct com- 
parison with dislocations in idealized random alloys. This 
is achieved by performing simulations for dislocations in 
Sio.sGeo.s alloys artificially constrained to have a random 
occupation of the lattice sites. We find that the disloca- 
tion formation energy (3.7 eV/A) is 2.0 eV/A greater 
than in the segregated alloy (1.7 eV/A). Surprisingly, 
however, the ideal random alloy has a smaller core ra- 
dius of 3.9 A with a correspondingly lower core energy 
of 0.31 eV/A, as determined by a linear fit of Eq. ([4]). 
The parameter fc of Eq. (|4]) is found to be 5.4 GPa. This 
value is, larger than the value found for the segregated 
alloy, as expected. 

To understand this we note that preventing atoms to 
exchange identities and occupy their energetically favored 
sites around the dislocation core results in weaker bonds 
in this area which are easier to deform. This explains 
the lower core energies associated with perfectly random 
alloy systems. However, despite higher core energies, the 
significant segregation outside of the core radius Rc (Fig. 
12) lowers the overall dislocation formation energy for seg- 
regated systems by reducing the overall elastic energy 
EeiiR). 

In summary, by partitioning the dislocation formation 
energy into elastic and segregation components we have 
established a general framework for calculating the core 
properties of dislocations in segregated alloys. This ap- 
proach opens up new possibilities for the study of more 
complex processes such as 'alloy drag' on the moion of 
dislocations [3], the infiuence of local segregation (in- 
cludin g v acancies) on the Peierls barrier and kink mo- 
bility [l6|, as well as dislocation- induced segregation in 
confined systems such as quantum dots. 
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